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Abstract 


The  characteristics  of  puised  diode  lasers  and  linear  CCD  array  opticai  detectors  are  explored  in  the 
context  of  high  speed  precision  dimension  and  location  measurements  by  optical  diffraction  and 
interference.  It  is  suggested  that  pixei  based  metrology  may,  with  appropriate  signal  extraction  and 
processing  methods,  be  made  mom  accurate  by  making  edges  /ess  weli  defined:  if  the  mechanism  of 
definition  ioss  is  weli  understood,  e.g.,  diffraction,  tfie  intensity  distribution  over  many  pixels  can  be 
matched  to  a  template  to  within  a  small  fraction  of  a  pixel,  obviating  the  need  to  make  sub-pixel 
interpolations  based  on  local  intensities.  The  specific  example  of  making  remote  flash  measurements 
of  the  diameter  and  location  of  an  opaque  (^iindrical  object,  e.g.,  a  fast-moving  textile  strand,  is 
examined.  A  model  is  developed  for  calculating  the  diffraction  pattern  given  the  target  and  instrument 
parameters.  The  results  are  examined  and  compared  with  experimental  data.  Algorithms  for  finding 
the  diameter  and  location  given  the  diffraction  and  interference  pattern  and  the  instrument  parameters 
are  described.  The  performance  statistics  of  a  prototype  instrument  are  measured  in  a  data  run  of 
1000  laser  shots.  Location  and  diameter  measurement  standard  deviations  are  0.13  pixel-widths 
(1.69  lun)  and  0.36  pixel-widths  (4.68  pm)  respectively,  implying  that  sub-micron  precisions  are 
obtainable  by  averaging  a  few  dozen  measurements.  Sources  of  drift  and  scatter,  and  their 
reduction,  are  discussed  in  the  context  of  designing  and  building  an  improved  second  generation 
instrument. 
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Introduction 


Pulsed  semiconductor  diode  lasers  and  CCD  detectors  invite  low  cost,  high  speed,  high  precision 
optical  diffraction  and  interference  measurements  of  mechanical  locations  and  dimensions.  A  diode 
laser,  a  high  current  pulsed  power  supply,  and  a  triggering  circuit  [8]  can  be  built  into  a  penlight-size 
module  capable  of  delivering  fast  (-  200  nsec)  pulses  of  high  power  (-  10  watts)  in  the  near  infra-red 
(«  0.9  pm)  at  high  repetition  rates  (»  2  kHz)  inexpensively  (-  $200).  A  linear  CCD  array  of  256 
elements,  its  docking  drcuits,  and  output  signal  amplifiers  can  be  purchased  as  an  evaluation 
board  [4]  or  user-assembled  into  a  compact  module  for  another «  $200.  The  dimensional  accuracy  of 
the  CCD  pixelation  (typically  on  13  pm  centers)  makes  it  an  attractive  ruler.  Given  these  inexpensive 
sources  and  detectors  the  designer  can  imagine  dozens  of  optical  measurement  configurations,  each 
adapted  to  the  needs  of  a  spedfic  application. 

The  dimensional  accuracy  of  CCD  pixelation  must  be  extremely  high:  the  only  practical  way  to 
maintain  layer-to-layer  alignment  of  large  scale  integrated  drcuit  devices  is  for  each  process  step  to 
be  dimensionally  accurate  in  an  absoiute  sense  [6].  On  the  other  hand,  despite  the  accuracy  of  the 
ruler,  in  the  absence  of  a  good  sub-pixel  interpolation  scheme  the  predsion  of  location  and  dimension 
measurements  made  with  a  13  pm  ruler  can  be  at  best  mediocre  relative  to  the  sub-pm  measurement 
requirements  of  pradical  interest  in  many  contexts.  Sub-pixel  interpolation  is  usually  done  by  a  local 
gray-level  method,  e.g..  a  gray-registering  pixel  straddled  by  a  white-registering  pixel  and  a  black¬ 
registering  pixel  is  defined  to  be  partially  obscured  by  a  step-function  edge  whose  location  is  linearly 
interpolated  via  the  gray-level  intensity  [9]. 

The  rrxxlels  and  experiments  described  in  this  report  suggest  that,  with  several  provisos,  CCD  based 
location  and  dimension  measurements  can  be  made  more  accurate  by  making  edges  less  well 
defined: 

•  the  method  of  edge  definition  loss  is  well  understood,  e.g.,  diffraction; 

•  the  parameters  that  dictate  the  dimensionai  details  of  the  spread-out  edge,  e.g.,  optical 
wavelength  and  apparatus  geometry,  are  known  with  adequate  accuracy; 

•  model-based  signal  extraction  and  processing  methods  are  utilized. 

Under  these  drcumstances,  the  intensity  distribution  of  an  edge  spread  out  by  diffraction  over  100  or 
more  pixels  can  be  matched  to  a  template  to  within  a  small  fraction  of  a  pixel.  I  suggest  that  the 
result  can  be  substantially  more  precise  than  sub-pixel  interpolation  based  on  iocal  gray-levels. 

In  this  report  I  describe  a  system  that  uses  a  pulsed  diode  laser,  a  CCD  evaluation  board,  and  a  few 
inexpensive  optical  components  to  make  flash  measurements  of  the  diameter  and  location  of  a 
stranded  material  composed  of  several  hundred  fibers  (•>  10  pm  fiber  diameter,  «  200  pm  strand 
diameter)  moving  at  high  speed  (>  50  m-secr>).  The  laser  and  the  CCD  are  controlled  by  a  single¬ 
board  microcomputer^  that  communicates  with  a  computer  workstation  that  runs  the  inversion 
algorithms  and  reports  the  results  and  statistics. 

The  harsh  environment  of  the  application  practically  precludes  an  imaging  system.  This  environment, 
and  a  requirement  for  unobstructed  operator  access  to  the  strand,  dictate  an  instrument  with  a  large 
standoff  (>  150  mm).  At  unit  magnification  with  this  standoff  an  image  would  be  only  about  15  pixels 
across.  Thus  even  carefully  implemented  gray-level  based  sub-pixel  interpolation  methods  would 


'Actually  ttw  apfMiatus  uses  several  lasers  and  CCDa  to  make  a  set  of  syrrchronous  dameter  and  location  measurements 
that  are  fused  to  give  a  measurement  of  higher  level  inisrest  e.g.,  for  process  diagnostics  and  control.  The  details  win  be 
discuased  in  another  report 


yield  relatively  low  precision.  Higher  than  unit  magnification  is  in  principle  a  solution,  but  it  would 
require  an  impractically  long  optical  system,  e.g.,  magnification  of  5x  with  a  150  mm  standoff  would 
put  the  detector  900  mm  (3  feet)  from  the  stramd.  Folded  optical  paths  are  of  course  a  possibility,  but 
prohibitively  expensive  in  the  application's  context 

Laser  diffraction  methods  for  measuring  diameters  have  developed  rapidly  in  the  last  few  years, 
primarily  in  response  to  the  process  control  requirements  of  the  fast  growing  optical  fiber 
communication  industry,  and  possibly  also  in  response  to  the  similar  requirements  of  specialty  metal 
wire  drawing  industries,  a.g.,  tor  integrated  circuit  wire  bonding.  The  optical  fiber  case,  at  least,  is 
dearly  amenable  to  diffraction  measurements  at  large  scattering  angles,  where  due  to  refraction  and 
internal  reflection  a  large  angle  signal  is  available  in  transmission  [2],  [3]. 

However  strands  composed  of  multiple  fibers  are  practically  opaque,  and  only  diffusely  reflecting, 
preduding  large  angle  scattering  methods.  Having  thus  ruled  out  imaging  and  large  angle  scattering, 
what  is  left  is  small  angle  scattering  in  the  forward  diredion,  i.o.,  analysis  of  the  shadow  that  the 
strand  casts  when  interposed  between  the  laser  and  the  detector.  This  shadow,  actually  a  Fresnel 
diffraction  pattern^,  contains,  I  will  show,  two  approximately  independent  sources  of  predse  (albeit 
computationally  expensive  to  extract)  information  about  the  strand  diameter  and  location.  The 
method  requires  virtually  no  optics:  an  inexpensive  (-  $40)  and  uncritical  microscope  objective 
simplifies  the  analysis  by  collimating  the  laser,  and  an  equally  inexpensive  and  uncritical  cylindrical 
lens  increases  signal-to-noise  by  focusing  the  longitudinal  direction  onto  the  linear  array.  This 
simpiidty  and  robustness  is  essential  in  the  application’s  environment,  which  indudes  severe 
vibration,  severe  electrical  interference  from  the  switching  of  kiloampere  currents,  high  temperature, 
humidity  in  excess  of  1 00%,  and  a  high  density  of  particulates  that  ding  tenadously  to  lenses  and 
other  glass  optical  elements.  These  constraints  preclude  a  finicky  optical  system:  the  predsion  of  my 
measurement  has  to  come  from  someplace  other  than  the  predsion  and  stability  of  the  mechanical 
and  optical  alignment. 


Theory  Behind  the  Model 

Modeling  the  shadow  of  an  opaque  wire-like  object  is  straightforward,  although  a  straightforward 
implementation  of  the  spatial  integrals  is  computationally  expensive.  Since  I  work  in  an  environment 
in  which  computational  expense  only  indirectly  translates  into  monetary  expense,  my  model  does  the 
integrals  by  brute-force.  This  slow-but-sure  approach  achieves  flexibility  and  an  absence  of  pitfalls 
that  might  not  be  achieved  by  a  clever  implementation,  e.g.,  one  employing  approximations  to  the 
phase  factors  to  facilitate  symbolic  rather  than  numerical  integration. 

Since  the  light  scatterer  is  wire-like,  I  will  model  the  problem  tor  cylindrical  waves,  so  that  all  planes 
normal  to  the  scatterer  are  equivalent,  and  only  one  needs  to  be  considered.  The  resulting  geometry 
of  the  scattering  (or  "source”)  plane  and  the  detector  plane  is  illustrated  in  Figure  1 .  A  coordinate  axis 
labeled  coinddes  with  the  wire-like  scatterer.  An  orthogonal  coordinate  axis  passes  through 
^trc  -  0-  opf'cal  electric-field  in  the  >,reZ,„:-plane  has  amplitude  and  phase  given  by  E(y^)  and 

in 


diffraction  patterns  are  broadly  daseed  ae  Fresnel  or  Fraunhoffer  based  on  Ihe  simpiifytng  approximations  that  can  be 
made  about  the  phase  factors  under  the  integral  signs  in,  e.g.,  Huygens’  Principle  (10)  or  Green's  Function  [7]  models  of  optical 
field  propagation.  For  the  purposes  of  this  report,  it  is  reasonable  to  say  that  a  Fresnel  pattern  is  a  shadow  (no  lens  used)  and 
a  Fraunhoffer  pattern  is  an  image  (ler«  used).  In  both  cases,  regions  that  in  geometrical  optics  would  contain  abrupt 
transitions  between  light  and  dark  acquire  oscillatory  structure  via  diffraction  and  interference  effects. 


3 


where  (d  is  the  optical  frequency.  The  wavenumber  t  =  S = ^,  where  c  is  the  speed  of  light  and  k  the 

Of^icai  wavelength,  is  the  magnitude  of  the  propagation  vector,  whk^  lies  along  a  third  axis  x 
orthogonal  to  the  yjrc^sre'P^^^  passing  through  its  origin.  The  polarization  vector  does  not  have 
to  be  specified,  but  to  be  conaete  let  us  say  it  is  in  the  direction.  A  detector,  realized  as  a  linear 
CCD  array,  lies  at  x  =  d,  along  coordinate  axis  which  is  parallel  to  y^.  At  any  point  on  the  y^ 
axis  the  optical  electric-field  is  given  by  a  Huygens’  Principle  integration  over  y^: 


^(ysrc) 
J* 


.[!+. 


^^*(ysrc-ydJ^ 

where  y^  and  y^  are  the  y^  coordinates  of  the  left  and  right  edges  of  the  strand. 
The  factors  of  this  equation  are  understood  as  follows: 


(1) 


1 .  ■  — :  The  electric-field  amplitude  at  attributable  to  the  Huygens’  wavelet 

originating  at  y^.  The  radiated  electric-field  falls  off  linearly  with  distance  between 
source  and  detector,  not  quadratically:  it  is  the  energy  density,  proportional  to  the 
square  of  the  electric-field,  that  falls  off  as  the  inverse  square  of  distance. 

2.  [  1  :  Kirchoff’s  modification  of  Huygens’  Principle,  the  "obliquity 

factor,”  which  assures  that  Huygens’  wavelets  propagate  only  in  the  forward  direction. 
The  factor  is  usually  written  l-t-oos8,  where  0  is  the  angle  of  wavelet  propagation  with 
respect  to  the  forward  direction. 

3  ^ * (♦  (ysre^ + *<ysrc~y<k/~““^  .  optical  phase  at  y^  lags  the  phase  at  y^  by  the 
wavenumber  times  the  path  distance  between  y^  and 


Photodetectors  in  general,  and  CCDs  specifically,  respond  to  incident  optical  energy,  which  is 
proportional  to  |E(y^)p  integrated  over  the  detector  area  and  multiplied  by  the  exposure  time.  The 

KB  C 

proportionality  constant  is  where  k  is  the  dielectric  constant  of  the  medium,  e.g.,  air,  and  is 

the  permittivity  of  free  space,  times  factors  describing  the  detection  efficiency  and  the  proportionality 
between  incident  energy  and  output  signal  voltage.  Since  all  these  factors,  and  also  the  pixel  height, 
will  ultimately  be  normalized  away,  I  will  simply  dispense  with  them  now  and  write,  for  the  signal  from 
the  A-th  pixel; 


««)  =  . 


M 

r: 

•V. 


m 


W. 


pixel 


where  y^,  is  the  y^,  coordinate  of  the  "left  edge"  of  the  n-th  pixel,  and  is  the  width  of  the 
unmasked  active  area  of  the  pixel.  Typical  values  are  13 pm  pixel-to-pixel  spacing,  and  5\m  active 
width. 


Implementation  of  the  Model 


Doing  the  integrals  over  y^  and  y^,  requires; 

•  replacing  the  integration  limits  of  -<»>  and  +<»  on  y^  with  appropriate  finite  limits; 
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ngura  1.  GEOMETOY  OF  THE  APPARATUS. 

Not  to  scale:  d  is  approximately  167  mm,  is  13  pm,  and  the  total  width  of  the  CCD  (256  calls) 

is  approximately  3.3  mm 

•  specifying  appropriate  forms  tor  E(y^  and 

•  specifying  appropriate  step  sizes. 

Many  of  the  choices  are  heuristic,  justified  by  the  absence,  when  they  are  made,  of  artifacts  in  the 
result,  robustness  of  the  result  against  small  changes  in  any  of  the  forms  and  parameters  selected, 
and  ultimately,  agreement  between  the  model  and  experimental  data.  The  following  short  paragraphs 
outline  a  successful  course. 


Doing  the  integral  over  requires  specifying  an  appropriate  step  size.  Some  computational 
efficiency,  especially  with  respect  to  being  able  to  reuse  partial  msults  from  example  to  example,  is 
greatly  enhanced  if  we  assume,  and  justify  by  comparison  between  the  resulting  calculation  and  the 
data,  that  the  optical  phase  changes  sufficiently  little  over  the  active  width  of  a  pixel  that  it  is  valid  to 
use; 


S(«)- 


W. 


pixd 


Using  this  approximation  is  equivalent  to  saying  that  the  phase  changes  slowly  enough  over  a  pixel 
width  that  it  is  adequate  to  average  the  electric-field  vector  components  over  a  pixel  and  take  the 
square-magnitude  of  the  result  instead  of  averaging  the  square-magnitude  over  a  pixel. 
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The  problem  of  'replacing  the  integration  limits  of  -«>  and  -t^o  on  with  appropriate  finite  limits”  is 
intimately  connected  with  "specifying  appropriate  forms  for  E(y^  and  "  The  steps  are: 

•  Assume  a  collimating  lens,  the  miaoscope  objective  mentioned  in  the  Introduction,  that 
transforms  the  laser  beam  pattern  into  plane  waves  in  the  scattering  plane:  thus 
=  consteua,  and  we  might  as  well  call  the  constant  zero; 


•  Recognize  that  ideal  laser  beams  have  a  gaussian  radial  structure,  so  we  would  like  to 
be  able  to  write: 


)■ 


where  is  a  scale  length  proportional  to  an  appropriate  transverse  linear  dimension  of 
the  active  laser  channel: 


•  Reason  that  the  appropriate  scaling  lengfo  is  given  by 


ln2 


tan  (sin  AW) 


where  R  is  the  half  aperture  (radius)  of  the  microscope  objective,  na  is  the  numerical 
aperture  of  the  microscope  objective,  and  'S  the  angular  full-width-at-half- 

maximum  of  the  laser  beam  divergence  perpendicular  to  the  long  dimension  of  the 
scatterer; 

•  Recognize  that  unless  the  integral  is  somehow  cut  off,  e.g.,  for  \y^^\  >  R,  the  integration 
is  both  physically  unreasonable  and  numerically  cumbersome.  Also,  discover  empirically 
that  unless  the  cutoff  is  made  smoothly,  artifacts  that  are  very  sensitive  to  the  details 
appear  in  the  result  Heuristically,  a  way  not  to  generate  artifacts  is  to  integrate  y^^ 
between  ~R  and  +R  while  smoothly  reducing  E(y^  to  zero  at  the  limits  by  parabolicaily 
truncating  the  gaussian,  i.e.,  taking: 

•  Rnally,  learn  empirically  that  step  sizes  for  the  integrals  over  y^  and  of  1  pm,  giving 
13  steps  per  inter-pixel  spacing,  and  5  steps  across  each  active  pixel  width,  are  a 
reasonable  compromise  between  smaller  step  sizes  resulting  in  excessive  computational 
time  and  larger  step  sizes  resulting  in  artifacts  such  as  spurious  modulation. 

Appendix  A  is  the  program  huygens  that  implements  this  model. 


Comparison  of  Model  with  Data 

In  this  section  I  compare  data  collected  with  prototype  apparatus  of  approximately  the  geometry 
described  in  Table  1  with  calculations  generated  by  the  huygens  program  using  these  parameters.  In 
the  following  figures  illustrating  the  model  the  calculation  is  degraded  to  6-bits  of  signal  amplitude 
resolution,  corresponding  to  the  6-bits  of  actual  resolution  available  from  the  analog-to-digital  flash 
converter  in  the  apparatus;  data  and  model  are  each  subjected  to  a  linearly  weighted  five  point  sliding 
average  before  plotting. 
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Tabtol 


Paramoters  of  the  apparatus  generating  the  data 
'  .4  of  the  mode}  simulating  the  data 


Parameter 

Value 

Comment 

CCD  cells 

256 

Number  of  cells  in  the  linear  CCD  detector 

Pixel  spacing 

13.0  pm 

Pixel-to-pixel  spacing 

Pixel  size 

5.0  pm 

Active  width  of  each  pixel 

Distance 

167  mm 

Target  to  detector  distance 

Wavelength 

0.895  pm 

Near  infra-red  pulsed  diode  laser 

Lens  radius 

2.5  mm 

Half-diameter  of  collimating  lens 

NA 

0.1 

Numerical  Aperture  of  collimating  lens 

FWHM 

30.0  deg 

Full-Width-at-Half-Maximum  laser  divergence 

ADC  range 

O-bits 

Resolution  of  analog-to-digital  converter 

Sources  of  error  that  should  be  borne  in  mind  when  comparing  data  and  model  include 

•  the  actual  wire  diameter  is  uncertain  because  on  the  one  hand  the  diameter  is  reduced 
by  stretching,  and  on  the  other  hand  it  is  increased  by  enamel  insulation; 

•  the  distance  is  uncertain  by  mechanical  measurement  error,  and  by  imperfect  collimation 
of  the  laser  beam,  which  may  be  slightly  convergent  or  divergent. 

Figure  2  compares  experimental  data  for  a  26  awg^  wire,  nominal  diameter  404  pm,  with  the  model 
for  a  target  diameter  of  400  pm.  This  diameter  is  about  twice  the  anticipated  strand  diameter  in  the 
application  of  interest.  The  graphs  for  data  and  model  are  slightly  offset  from  each  other  horizontally 
to  correct  for  the  wire  being  inexactly  centered  on  the  CCD.  The  contrast  in  the  data  is  somewhat 
less  than  the  contrast  in  the  model,  due  no  doubt  to  stray  light,  scattering  of  the  laser  beam  by  optical 
surfaces,  optical  imperfections,  dust,  etc.  The  correlation  between  data  and  model  are  clearly 
extremely  high. 

Figure  3  compares  experimental  data  for  a  28  awg  wire,  nominal  diameter  320  pm,  with  the  model  for 
a  target  diameter  of  320  pm.  This  diameter  is  about  50%  larger  than  the  anticipated  strand  diameter 
in  the  application  of  interest.  The  remarks  made  with  respect  to  Figure  2  again  apply.  Systematic 
diameter  related  differences  between  Figures  2  and  3  are  easily  discerned.  For  example,  the 
prominent  wiggle  between  each  toe  and  shoulder  of  the  shadow  has  moved  outward  in  both  data  and 
model.  But  the  wiggle  that  in  Figure  2  is  midway  between  each  toe  and  shoulder  of  the  shadow  has 
in  both  the  data  and  model  of  Figure  3  moved  all  the  way  up  to  each  shoulder  in  the  model,  but  only 
part  way  there  in  the  data.  We  can  probably  conclude  that  the  26  awg  wire  is  slightly  oversize,  i.e.,  in 
the  direction  of  26  awg.  Since  the  thickness  of  the  enamel  insulation  has  not  been  taken  into 
account,  this  is  reasonable. 

Figure  4  compares  experimental  data  for  a  32  awg  wire,  nominal  diameter  202  pm,  witti  the  model  for 
a  target  diameter  of  approximately  200  pm.  This  diameter  is  about  equal  to  the  anticipated  strand 


^American  Wire  Gauge  (awg)  dameterB  increase  (ogarithmically  with  decreasing  gauge  number.  In  this  report  the  tabular 
relatioftship  (1)  is  repressed  by  dituneur(itm)  •  ap(mawg  +  b)  where  m  >  (/oj (79.87/404.0))/ (40.0 -26.0)  and 
b  «  /oj  ^9.87)- 40.0  xm 
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Flgura  2.  26  AWG  WIRE  DATA 

Top,  data  for  26  awg  wira,  approximately  404  pm  diameter; 
bottom,  model  for  400  pm  diameter. 

diameter  in  the  application  of  interest.  As  in  Figures  2  and  3,  except  for  degraded  contrast  in  the 
experimental  data  the  agreement  between  data  and  model  is  excellent.  The  wiggle  that  moved 
outward  between  Figure  2  and  3  has  moved  still  farther  out,  somewhat  more  so  in  the  model  than  in 
the  data.  This  probably  indicates  that  the  32  awg  wire  is  also  somewhat  oversize,  again  probably  due 
to  the  enamel  insulation. 

Figure  5  compares  experimental  data  for  a  38  awg  wire,  nominal  diameter  101  pm,  with  the  nxidel  for 
a  target  diameter  of  approximately  100  pm.  This  diameter  is  about  half  the  anticipated  strand 
diameter  in  the  application  of  interest.  The  trends  identified  in  Rgures  2,  3,  and  4  are  seen  to 
continue  in  a  smooth  and  predictable  fashion. 

There  are  three  obvious  conclusions  of  this  comparison: 

1 .  the  model  and  the  data  agree  well; 

2.  the  data  are  rich  in  diameter  dependent  feabires; 

3.  over  the  range  of  target  diameters  between  100  and  400  pm  the  dominant  feature  of  the 
data,  i.e.,  the  central  dip  that  intuitively  might  be  associated  with  the  geometrical 
shadow,  in  fact  changes  only  slightly  in  width;  thus  the  width  of  this  dip  is  not  a 
straightforward  measu.e  of  target  diameter; 

4.  the  wiggle  seen  midway  between  each  toe  and  shoulder  of  the  dip  in  Figure  2,  diameter 
-  400  pm,  has  by  Figure  5,  diameter  -  100  pm,  moved,  quite  systematically,  halfway 
across  the  CCD. 
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Flgur«3.  28  AWG  WIRE  DATA 

Top,  data  for  28  awg  wire,  appioximataly  320  |un  diamatar; 
txMom,  modal  for  320  |tm  diamatar. 


Rgure  6  aids  in  understanding  the  two  distinct  diameter  dependent  features  in  the  shadow  by 
showing  the  model  result  for  a  200  pm  diameter  target  with  the  distance  parameter  set  to  1  mm,  10 
mm,  and  100  mm  respectively.  These  three  shadows  can  also  be  compared  with  Rgure  4,  which 
corresponds  to  the  prototype  apparatus  distance  of  167  mm.  In  Rgure  6,  as  in  Figures  2-5,  the 
model  has  been  subjected  to  the  same  6-bit  digitization  and  five-point  linearly  weighted  smoothing  as 
the  data.  At  1  mm,  the  top  of  Figure  6,  the  width  of  the  central  dip  is  essentially  the  diameter  of  the 
target^,  i.e.,  it  is  essentially  a  geometrical  shadow  with  some  fine  structure  detemnined  by  diffraction 
and  interference.  At  10  mm,  the  middle  of  Figure  6,  the  width  of  the  central  dip  has  increased 
significantly,  and  diffraction  effects  extend  for  several  tens  of  pixel-widths  to  each  side  of  each 
shoulder.  At  100  mm,  the  bottom  of  Rgure  6,  diffraction  and  interference  effects  extend  at  least  a 
hundred  pixel-widths  beyond  each  shoulder. 

It  is  especially  apparent  in  the  bottom  member  of  Rgure  6,  as  well  as  in  Figure  5,  that  there  are  two 
phenomena  contributing  to  the  structure  of  the  data.  Most  of  the  signal  power  is  contained  in  a 
pattern  of  wiggles  whose  spacing  and  amplitude  decrease  systematically  with  distance  from  the 
shoulders  of  the  geometrical  shadow.  The  peak  and  valley  locations  in  this  structure  do  not  obviously 
depend  on  the  wire  diameter,  although  there  is  a  distinct  dependence  of  the  oscillation  amplitude  on 
the  diameter.  Superimposed  on  this  weakly  diameter  dependent  structure  is  a  pattern  of  uniformly 
spaced  wiggles  whose  spacing  increases  dramatically  as  the  diameter  decreases.  This  second 
feature  contains  regrettably  little  signal  power,  and  thus  while  it  is  easily  visible  in  the  model,  parts  of 
it  are  sometimes  difficult  to  identify  in  the  noise  of  the  experimental  data. 


^The  measurement  can  be  verified  simply  by  multiplying  the  ratio  of  the  dip  widtti  to  the  figure  width  by  the  actual  width 
representod  by  the  figure  width,  2S6  CCD  channels  times  13  pm  per  channel. 
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FIgura  4.  32  AW6  WIRE  DATA 

Top.  data  for  32  awg  wire,  approximately  202  pm  diameter; 
bottom,  model  for  200  pm  diameter. 

The  oondusion  of  these  experiments  is  that  there  is  a  wealth  of  diameter-sensitive  information  in  the 
data.  The  question  now  is  whether  the  inverse  problem,  determining  the  diameter  from  the  data,  is 
possible,  and  if  possible,  if  it  is  practical.  Answering  these  questions  affirmatively  is  aided  by  an 
understanding  of  the  physical  origins  of  the  features  this  section  has  identified. 


Building  Intuition  about  the  Method 

The  real  and  simulated  data  are  essentially  iden^cal.  We  conclude  that  we  know,  at  least  to  the  level 
of  our  measurement  precision,  how  to  calculate  what  the  apparatus  will  tell  us.  I  will  now  show  that 
we  can  also  understand  the  physical  origin  of  the  major  features  of  the  data  as  characterized  when 
the  target-to-CCO  distance  exceeds  a  few  millimeters; 

1 .  There  is  a  dip  in  light  intensity  that  is 

•  centered  on  the  geometrical  shadow,  and 

•  much  broader  than  the  geometrical  shadow; 

2.  The  center  of  the  dip  is  somewhat  brighter  than  the  region  a  little  off  center; 

3.  The  light  intensity  on  each  side  of  the  dip  is  oscillatory  with  a  broad  structure  that: 

•  seems  to  the  eye  to  be  independent  of  strand  diameter,  (in  spacing,  although  not 
so  independent  in  intensity),  and 


FIgurt  5.  38  AWQ  WIRE  DATA 

Top,  data  for  awg  38  wira,  approximately  101  pm  diameter; 
bottom,  model  for  1 00  pm  diameter. 

•  in  which  the  spatiai  frequency  of  the  oscillations  increases  (the  spaong 
decreases)  with  distance  from  the  strand  center; 

4.  Superimposed  on  the  oscillatory  structure  is  a  fine  structure  that  is: 

•  extremely  sensitive  to  strand  diameter,  and 

•  whose  spatial  frequency  (and  thus  the  distance  between  manifestations)  is  a 
constant  that  depends  inversely  on  diameter. 

Of  these  features,  the  third  one  (the  broad  oscillations  of  increasing  spatial  frequency)  is  most 
important  for  the  interim  inversion  procedure  I  have  implemented.  The  fourth  feature,  the  fine 
structure  of  constant,  strongly  diameter  dependent  spacing,  promises  to  be  most  important  for  future, 
higher  precision  implementations. 

Insight  into  the  origin  of  the  broad  oscillations  is  obtained  by  studying  the  shadow  of  a  knife  edge,  the 
bottom  of  Figure  7.  The  top  of  Figure  7  is  the  same  as  the  model  result  previously  shown  in  Figure  5. 
The  broad  oscillations  of  increasing  spatial  frequency  seen  in  the  data  and  model  are  thereby  easily 
recognized  as  the  shadows  of  the  edges  of  the  wire-like  target:  the  primary  structure  is  due  to 
diffraction,  and  it  deperis  on  each  edge  essentially  independently.  The  horizontal  offset  between  the 
top  and  bottom  parts  ot  -igure  7  is  essentially  the  measure  of  the  right  edge  location  with  respect  to 
the  wire  center.  The  offset  shown  was  obtained  by  visually  aligning  the  top  and  bottom  parts  of  the 
figure.  The  offset  measures  2  mm  in  a  total  of  128  mm  figure  width,  corresponding  to  256  pixel- 
widths  of  13  pm  each,  or  a  wire  radius  of  52  pm,  compared  with  the  nominal  value  of  50  pm. 


ngur*  6.  DIAMETER  DEPENDENCE  MODEL 

The  model  for  200  (im  diameter  for  target-to-CCD  distance  of 
(op,  1  mm;  middle,  fO  mm;  and  bottom,  100  mm. 

The  remaining  fine  structure  is  an  interference  effect  that  depends  on  both  edges  simultaneously. 
Rgure  8  shows  the  model  result  for  two  1  pm  wide  strips  200  pm  apart,  i.e.,  a  two  slit  interference 
pattern  with  a  slit  separation  the  same  as  the  target  diameter,  compared  with  the  nxxlel  result  for  a 
200  pm  diameter  wire.  Aithough  generated  for  this  demonstration  by  the  huygens  program,  the 
two-slit  part  of  Figure  8  is  really  just  an  offset  cosine  function  whose  wavelength  is  the  product  of 
distance  and  optical  wavelength  divided  by  target  diameter.  Comparing  it  with  the  model  for  the  200 
pm  wire,  it  is  easily  seen  that  each  peak  in  the  two-slit  pattern  corresponds  precisely  to  an 
interference  "wiggle”  perturbing  the  broad  diffraction  structure  of  the  wire.  The  fine  structure, 
because  it  depends  on  both  edges  simultaneously,  is  exquisitely  sensitive  to  the  target  diameter.  But 
this  sensitivity  is  difficult  to  exploit  because  the  signal  power  associated  with  it  is  small  relative  to  the 
signal  power  associated  with  the  broad  oscillations. 
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Figurar.  DIFFRACTION  SHADOWS 

Diffraction  shadow  of  bottom,  a  knife  edge,  and  Ipp,  a  100  diameter  wire. 


The  intuitive  picture  I  have  given  is  piausible,  and  the  agreement  between  the  brute  force  numericai 
integrations  and  the  data  support  the  daimed  deep  understanding.  While  it  is  difficult  to  predict  the 
relative  intensities  of  the  diffraction  and  interference  effects  by  symbolic  calculation,  it  is  nevertheless 
straightforward  to  predict  the  spadng  of  the  two  kinds  of  osdilatory  terms.  I  condude  this  section  with 
simple  derivations  of  the  spadng  of  osdilations  assodated  with  diffraction  at  a  knife  edge,  and  with 
interference  between  two  such  diffracted  waves. 

The  situation  at  the  knife  edge  is  illustrated  by  Rgure  9.  A  plane  wave  is  moving  from  left  to  right.  It 
is  interrupted  by  the  knife-edge  toward  the  left  side  of  the  picture.  In  the  upper  half  of  the  piclire  the 
plane  wave  continues  on  to  the  detector  plane  at  the  far  right.  A  cylindrical  Huygens’  wavelet 
propagates  into  both  halves  of  the  space  to  the  r^ht  of  the  knife-edge.  The  results  are 

•  The  cylindrical  wavelet  canies  energy  into  the  region  that  would  be  in  the  dark  shadow 
space  were  optics  geometrical; 

•  The  interference  between  the  plane  waves  and  the  cylindrical  wave  produces  an 
osdilatory  structure  in  the  region  that  would  be  in  the  completely  unshaded  space  were 
optics  geometrical. 

The  geometry  is  simple:  constructive  interference,  resulting  in  increased  light  intensity,  occurs  in  the 
vidnity  of  regions  where  the  path  difference  between  plane  and  cylindrical  waves  is  an  even  number 
of  half  wavelengths;  destructive  interference,  resulting  in  decreased  light  intensity,  occurs  in  the 
vidnity  of  regions  where  the  path  difference  is  an  odd  number  of  half  wavelengths.  For  the  small 
angles  that  are  applicable  the  path  difference  between  plane  and  cylindrical  waves  reaching  is 


Figure  8.  TWO  SUT  INTERFERENCE  PATTERN 

Interference  bottom,  between  two  1  )im  wide  slits  200  |un 
apart,  compared  with  top  model  for  a  200  pm  diameter  wire. 


petti  dHtannoe 


Figure  9.  KNIFE-EDGE  GEOMETRY 
Geometry  for  calculating  knife-edge  diffraction  pattern. 


=  ^3- 

for  to  the  right  of  a  right  edge,  and 
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for  to  the  left  of  a  left  edge,  where  d  is,  as  previously,  the  distance  between  the  target  plane  and 
the  detector  plane.  To  the  left  of  a  right  edge  and  to  the  right  of  a  left  edge  there  is  only  the  cylindrical 
wave,  no  plane  wave,  so  there  are  no  interference  osdllations  in  these  regions,  only  a  rapid  falloff 
with  angle  corresponding  more-or-less  to  the  l-cos  d  obliquity  factor. 


The  maxima  and  minima  of  the  oscillatory  structure  outside  the  geometrically  shaded  region  occur 
whenAp  =  aV2,  at 

=  V2dnX  (2) 

where  Ay^  symbolizes  either  or  (y^-y^,)  as  the  case  may  be.  Even  values  of  n  are  the 

peaks,  odd  values  of  n  are  the  valleys,  measured  from  the  location  of  the  geometrical  shadow  in  the 
detector  plane.  Measuring,  more  naturally,  from  the  center  of  the  shadow  of  a  wire-like  target,  half 
the  target  diameter  is  added  to  A>^  for  locations  outside  the  geometrical  shadow,  and  Ay^  is 
subtracted  from  half  the  target  diameter  fbr  locations  inside  the  geometrical  shadow.  Measurements 
of  the  peak  and  valley  locations  on  data  and  simulations  agree  well  with  this  model. 


This  treatment  re-affirms  that  the  primary  osciilatory  structure  is  a  weak  source  of  target  diameter 
information,  in  that  it  is  only  the  offset  of  this  structure  from  target  center  location  that  measures 
diameter.  On  the  other  hand,  it  is  a  better  measure  than  would  be  the  location  of  a  geometrical 
shadow,  in  that  the  edge  iocation  information  is  spread  over  many  pixels,  and  can  be  measured  to 
within  a  small  fraction  of  a  pixel-width  by  a  mask  or  template  matching  algorithm.  In  contrast,  the 
iocation  of  a  geometrical  shadow  with  sub-pixel-width  precision  has  to  be  inferred  from  the  gray  level 
interpolation,  a  method  that  is  probably  in  general  less  precise. 


The  geometry  of  two  slit  interference  is  illustrated  by  Figure  10.  We  consider  the  interference  of  the 
two  cylindrical  wavelets  originating  at  the  left  and  right  edges  of  the  target.  This  is  the  Young’s  two 
slit  interference  pattern  problem  for  slit  widths  much  smaller  than  slit  spaa'ng,  where  the  *slit  spadng” 
is  our  target  diameter.  For  distances  Ay^,  not  too  far  from  the  center  of  the  interference  pattern 

d  D  ~  2D  ^ 

where  in  this  case  the  target  diameter  D  is  an  explicit  factor  in  the  result.  Peaks  and  valleys  are 
located  at  integer  values  of  n,  peaks  at  even  n  and  vaileys  at  odd  n.  What  we  see  in  reality,  when  this 
interference  pattern  is  superimposed  on  the  edge  diffraction  patterns,  is  a  wiggle  at  each  even  value 
of  n.  At  n  =  0,  i.e.,  at  the  center  of  the  diffr^tion  pattern,  interference  is  also  constructive;  this 
accounts  for  the  bright  spot  at  the  center  of  e^  pattern,  analogous  to  the  well  known  bright  spot 
found  at  the  center  of  the  shadow  of  a  small  disk  or  sphere  [5].  The  locations  of  the  regularly  spaced 
small  wiggles  in  the  real  and  simulated  data  are  in  excellent  agreement  with  the  prediction  of 
Equation  3  fbr  consecutive  even  values  of  n. 


This  interference  structure  is  a  tantalizing  source  of  diameter  measuring  information  for  at  least  two 
reasons; 

1 .  The  peaks  and  valleys  have  a  fixed  spacing,  corresponding  to  a  fixed  spatial  frequency; 
their  fourier  transform  is  a  uniformly  spaced  line  sp^rum,  in  contrast  to  the  continuum 
spectnjm  generated  by  the  edge  diffraction  pattern; 

2.  The  spadng  depends  inversely  on  the  diameter  D,  with  correspondingly  increasing 
relative  sensitivity  to  smailer  diameters. 

The  prototype  instrument  was  built  with  neither  sufficient  shot-to-shot  laser  pattern  reprodudbility  nor 
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RguralO.  TWO  SUT  GEOMETRY 
G«omatry  for  calculating  two  slit  intaiferance  pattern. 

sufficient  anaiog-to-digitai  conversion  resolution  to  be  able  to  demonstrate  this  potential  for  precision 
interferometric  measurement  of  diameter.  I  hope  to  remedy  this  in  the  future. 

Some  final  words  of  caution  are  in  order.  Rrst,  the  treatments  of  both  single  edge  diffraction  and  two 
edge  interference  ignore  issues  of  phase  shifts  on  scattering;  these  complicated  effects,  which 
depend  on  the  materials  involved,  e.g.,  may  be  different  for  metals  and  insulators,  sometimes  have 
effects  such  as  reversing  the  roles  of  even  and  odd  values  of  n  in  Equations  2  and  3.  Second,  strictly 
speaking  it  is  incorrect  to  separate  single  edge  diffraction  and  two  edge  interference;  there  is  a  single 
scattering  process,  more-or-less  described  by  Equation  1 .  Because  the  detection  process  destroys 
phase  information,  it  is  fundamentally  impossible  to  dissect  the  result  into  two  or  more  separate 
effects.  Nevertheless,  if  one  does  not  ask  too  many  questions  about  details,  particularly  details  of  the 
intensity  distribution,  but  instead  restricts  the  discussion  to  the  location  of  significant  features  such  as 
peaks  and  valleys,  this  treatment  is  heuristically  valid.  In  practice,  the  intuition  generated  by  the 
model  that  separates  single  edge  diffraction  and  two  edge  interference  was  essential  to  disoovering 
an  algorithm  to  extract  the  diameter  and  location  information  from  the  data. 


The  Diameter-Center  Algorithm 

The  diameter  and  center  are  obtained  via  an  algorithm  that  independently  locates  the  left  and  right 
edges.  The  diameter  is  defined  as  the  difference  in  edge  locations  and  the  center  is  defined  as  the 
mean  of  edge  locations.  Tbe  algorithm  operates  in  several  steps; 


1.  TTie  digitized  CCD  data  are  smoothed,  arnj  in  the  process  converted  from  integer^  to 
floating  point  values,  by  a  five  point  linearly  weighted  sliding  average,  as  illustrated  in 
Rgure  1 1  tor  a  30  awg  wire  target. 


Figure  11.  30  AWQ  VWRE  DIGiTIZED  DATA 

Top,  raw  digitized  CCD  data,  30  awg  wire,  approximately  254  pm  diameter 
(transmitted  as  twos-oompiement,  therefore  down  is  light,  up  is  dark). 
Bottom,  same  data  after  five  point  linearly  weighted  sliding  average. 


the  present  64>it  implementalion,  and  the  anticipated  S-txt  implementation,  it  is  efficient  to  actually  store  the  raw  data  as 
character  (byte)  sized  variables. 
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2.  The  smoothed  data  are  subtracted  from  a  stored  floating  point  background^  or  no  wire 
pattern  obtained  by  averaging  the  pattern  obtained  from  several  («  16)  laser  shots  as 
illustrated  in  Figure  12  which  shows  the  stored  background,  the  data  after  background 
subtraction,  and  the  model  result  for  a  260  pm  wire. 


Figure  1 2.  30  AWG  WIRE  PROCESSED  DATA 

Top,  stored  background  {no  wire)  file. 

Center,  data  of  Figure  1 1  after  subtraction  from  the  background. 
Bottom,  model  result  for  260  ^m  diameter  wire. 


*The  data  are  subtracted  from  the  background  rather  than  vice  versa  because  the  digitizer  is  set  up  to  deliver  twos- 
oomplement  data. 
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3.  The  background-subtracted  data  are  convolved  with  a  mask  or  template  that  is  a 
heuristically  selected  portion  of  the  diffraction  shadow  of  a  knife  edge,  Figure  7;  the 
portion  that  is  used  begins  two  pixels  to  the  right  (higher  pixel-numbers)  of  center  and 
extends  for  64  pixels.  These  convolutions,  and  the  convolutions  with  the  mirror-image  of 
the  mask,  are  illustrated  in  Rgure  13. 


Rgura  13.  30  AWG  WIRE  CONVOLUTIONS 

Top,  convolutions  of  a  portion  of  the  mask  of  Figura  7 
with  tha  smoothad,  normalizad  data  of  Figura  12. 

Bottom,  sama,  with  a  miiTor-imaga  of  tha  mask. 

Tha  heavy  vertical  bars  mark  nominal  edge  locations  found  by  tha  algorithm. 

4.  The  absolute  minimum  of  the  convolution  is  found,  and  then  the  first  maximum  to  the 
right  of  the  minimum  is  found;  this  maximum  is  marked  by  a  heavy  vertical  bar  in  Figure 
13. 

5.  A  polynomial  is  fit  to  a  seven-pixel-wide  region  centered  on  this  maximum. 

6.  The  right  edge  is  defined  as  the  fractional  pixel-number  at  which  this  polynomial 
maximizes. 

7.  The  left  edge  is  found  by  a  procedure  tfiat  is  the  mirror-image  of  the  four  previous 
steps,  as  shown  in  the  bottom  part  of  Figure  13. 

8.  The  diameter  and  center  are  obtained  by  respectively  subtracting  and  averaging  the 
two  edge  locations. 

This  procedure  fails  if  the  right  edge  falls  within  64  pixels  of  the  right  end  of  the  CCD  or  the  left  edge 
falls  within  64  pixels  of  the  left  end  of  the  CCD;  in  either  of  these  cases  an  alternative,  less  accurate, 
procedure  is  invoked.  The  alternative  procedure  fits  a  polynomial  to  a  seven-pixel-wide  region 
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centered  on  the  minimum  of  the  convolutions.  In  a  nx>derate  range  of  diameters  this  minimum  has  an 
approximately  fixed  ofteet  from  the  maximum  actually  sought 

Even  the  alternative  procedure  fails  if  an  edge  is  off  the  CCD;  when  this  happens  the  particular 
measurement  is  abandoned. 

At  the  end  of  the  procedure  a  preliminary  calibration  correction  is  invoked  by  a  procedure  that  adjusts 
the  diameter  via  a  low  order  polynomial  mapping  of  the  diameter-center  plane.  The  polynomial 
coefficients  are  found  by  running  the  algorithm  on  simulated  data  generated  by  the  huygens  program. 
The  largest  effect  of  the  correction  is  to  subtract  a  constant,  approximately  12  pixel-widths. 

Appendix  C  is  the  function  diacens  (and  several  functions  that  it  uses)  that  implements  this  model, 
and  Appendix  B  is  the  Indude  file  of  instrument  constants,  program  parameters,  and  structure 
definitions  needed  to  understand  parts  of  Appendix  C. 

We  end  this  section  by  explicitly  noting  that  the  background  subtraction  procedure  is  heuristic;  on 
strict  theoretical  grounds  it  is  wrong.  Equation  1  gives  the  electric-field  amplitude  and  phase  in  the 
detector  plane  as  an  integral  of  the  electric-rteld  amplitude  and  phase  over  the  un-ocduded  portion  of 
the  target  plane.  The  background  or  no  wire  picture  is  essentially  the  square  magnitude  of  this 
integral  when  there  is  no  occlusion  of  any  portion  of  the  target  plane.  But  taking  the  square 
magnitude,  i.e.,  translating  from  vector  arnplitudes  (which  add)  to  scalar  intensities  (which,  for 
coherent  illumination,  cannot  properly  be  added,  but  which  are  what  the  detector  reports) 
permanently  destroys  the  phase  information  that  is  needed  to  do  a  proper  background  correction.  In 
fact,  if  the  detection  process  did  not  destroy  this  phase  information,  the  integral  of  Equation  1 ,  with 
subscripts  det  and  src  interchanged,  could  be  used  to  do  the  inversion  exactly.  The  result  would  be 
the  laser  intensity  pattern  in  the  un-ocduded  part  of  the  target  plane,  and  blackness  across  the 
diameter  of  the  target.  As  a  practical  matter,  the  procedure  I  have  described  is  effective  despite  toe 
fact  that  it  casually  treats  intensities  as  if  they  could  engage  in  toe  arithmetic  that  nature  in  fact  has 
reserved  for  amplitudes. 


Instrument  Performance 

The  raw  data,  their  trend  line,  and  the  standard  deviation  of  their  distribution  for  1000  diameter-center 
measurement  made  in  a  single  run  lasting  between  10  and  1 1  hours  are  shown  in  Figures  14  and  15. 
The  target  was  a  30  awg  wire. 

The  mean  of  the  raw  diameter  measurements  is  21 .504  pixel-widths,  or  279.5  pm,  in  contrast  to  the 
actual  diameter  of  254.2  pm;  a  multiplicative  calibration  factor  having  no  significant  bearing  on  any 
performance  measures  would  correct  this  disaepancy.  The  spread  in  toe  raw  measurements  is  from 
a  maximum  of  22.67  pixel-widths  to  a  minimum  of  20.66  pixel-widths;  most  of  the  skew  in  the 
distribution  is  due  to  the  trend,  which  could  (but  will  not)  be  subtracted  out.  The  standard  deviation  is 
0.36  pixel-widths,  which  is  the  appropriate  error  to  assign  to  a  single  measurement.  The  relative 
standard  deviation  is  1 .67%.  Averaging  n  measurements  would  reduce  these  values  by  Vn,  e.g.,  toe 
average  of  100  measurements  would  be  expected  to  be  correct  within  0.036  pixel-widths,  or  0.167%. 

The  mean  of  toe  raw  center  measurements  is  pixel-number  128.00;  since  there  is  no  independent 
measure  of  toe  accuracy  of  toe  centering  this  value  is  entirely  nominal;  its  only  interest  is  as  a  fiducial 
value  with  respect  to  which  trend  and  scatter  can  be  evaluated.  The  spread  in  toe  raw 
measurements  is  from  pixel-number  127.57  to  pixel-number  128.38;  as  in  the  case  of  toe  diameter 
measurement,  most  of  toe  skew  in  toe  distribution  is  due  to  the  trend,  which  could  (but  will  not)  be 
subtracted  out.  The  standard  deviation  is  0.1 31  pixel-widths.  This  is  toe  appropriate  error  to  assign 
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Hgur«14.  DIAMETRIC  SCATTER  PLOT 
Scattor  plot  for  1 000  diamoter  measurements  over  approx.  1 0  hrs. 

to  a  Single  measurement  Averaging  n  measurements  would  reduce  this  value  by  Vn,  e.g.,  the 
average  of  1 00  measurements  would  be  expected  to  be  correct  within  0.0131  pixel-widths. 


Figure  IS.  CENTRIC  SCATTER  PLOT 
Scatter  plot  for  1000  center  measurements  over  approx.  10  hrs. 

For  both  the  diameter  and  the  center  measurements  the  trend,  or  long  term  drift,  accounts  for  about 
20%  of  the  the  total  error,  and  the  remaining  80%  is,  to  visual  inspection,  reasonably  random  shot-to- 
shot  noise.  Both  the  long  term  drift  and  the  shot-to-shot  noise,  I  believe,  are  due  primanly  to  the 
laser.  Drift  in  the  mean  laser  pattern  (the  stored  background  or  "no  wire"  pattern,  Rgure  12),  which  in 
the  present  implementation  is  recorded  only  once  before  the  start  of  data  collection,  I  believe 
accounts  for  most  of  the  trend.  Shot-to-shot  fluctuation  in  the  laser  pattern  is  substantial,  as  much  as 
25%  in  some  parts  of  the  pattern,  presumably  because  the  laser  is  operated  fairly  close  to  threshold, 
and  I  believe  that  the  shot-to-shot  noise  is  almost  entirely  due  to  this  fluctuation. 

The  long  term  drift  could  be  substantially  reduced  by  periodically  moving  the  wire  out  of  the  laser 
beam  and  revising  the  stored  background  pattern;  this  could  be  automated  easily.  Alternatively,  if  the 
mean  value  of  the  measurement  were  known  via  other  measurements  to  be  time  invariant,  the  trend 
line  could  simply  be  subtracted  out.  The  background  revision  alternative  is  obviously  preferable. 

The  shot-to-shot  noise  coukJ  be  substantially  reduced  by  introducing  a  spatial  filter  (focusing  lens, 
pinhole,  and  re-collimating  lens);  the  pattern  emerging  from  the  spatial  filter  would  continue  to  show 
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globai  intensity  fluctuations,  but  its  spatial  pattern  would  be  stabilized  to  essentially  the  gaussian 
TEMqq  mode.  Furthermore,  to  compensate  for  the  light  lost  in  the  spatial  filter  the  laser  would  have  to 
be  operated  at  higher  power,  thus  redudng  the  globai  intensity  fluctuations. 

The  noise  is  so  strongly  dominated  by  these  laser-related  factors  that  it  is  impossible  to  say  what  are 
the  next  largest  sources  of  technical  noise.  It  seems  likely  that  adding  a  spatial  filter  will  reduce  the 
drift  and  noise  levels  by  an  order-of-magnitude  or  more. 


Summary  and  Next  Steps 

The  diffraction  shadow  examined  in  this  report  can  be  productively  though  of  as  containing  two 
distinct  sources  of  metrological  information:  the  independent  shadows  of  the  right  and  left  edges,  and 
the  interference  pattern  between  the  two  diffraction  patterns.  The  shadows  have  been  shown  to  be 
suitable  for  finding  the  edges  via  a  mask  or  template  matching  operation  with  the  calculated  (or 
alternatively,  measured)  shadow  of  a  knife  edge.  This  method,  since  it  uses  the  information  from  the 
two  edges  independently,  is  error  prone  and  of  relatively  low  precision  compared  with  what  is 
potentially  available  in  the  interference  pattern,  whose  spatial  frequency  depends  explicitly  and 
delicataly  on  the  diameter  itself.  However  substantially  more  signal  power  is  available  in  the  edge 
shadows  than  in  the  interference  pattern.  With  tiie  low  signal-to-noise  and  coarse  analog-to-digital 
converter  implemented  in  the  prototype  described,  using  the  edge  shadows  has  proven  to  be  the 
expedient  alternative. 

Using  the  edge  shadow  method,  the  prototype  apparatus  has  demonstrated  the  ability  to  make 
diameter  measurements  with  a  standard  deviation  (expected  error  of  a  single  measurement)  of  0.36 
pixel-widths,  or  4.7  pm,  for  wire-like  targets  in  the  100  to  400  |im  diameter  range.  The  corresponding 
standard  deviation  for  the  center  location  measurement  is  0.13  pixel-widths,  or  1.7  pm.  The  laser 
flash  duration  is  approximately  200  nsec,  the  CCD  readout  time  could  be  as  short  as  100  nsec  per 
channel,  or  25.6  psec  per  line;  video  flash  digitizers  that  could  follow  this  10  MHz  data  rate  are  easily 
available.  Thus,  computation  times  aside,  measurements  could  be  repeated  at  a  rate  approaching 
40,000  sec*.  While  a  general  purpose  computer  could  not  begin  to  keep  up  with  this  repetition  rate, 
spedai  purpose  hardware,  e.g.,  an  array  processor  card,  could  easily  come  within  an  order  of 
magnitude  of  keeping  up  in  real  time.  When  n  measurements  are  averaged  the  uncertainty  of  the 
average  is  the  uncertainty  of  a  single  measurement  divided  by  Vn;  thus  with  a  0.25  second  averaging 
time  we  could  ideally  obtain  a  factor  of  1000  improvement,  /.e.,  4.7  nm  diameter  measurement  and 
1 .7  nm  center  location  measurement.  In  short,  sub-pm  precision  is  easily  obtainable  by  averaging  a 
few  dozen  measurements. 

A  more  interesting  alternative  to  the  brute-force  approach  of  making  and  averaging  many 
measurements  in  a  short  time  is  to  look  for  opportunities  to  extract  and  use  the  interference 
information  that  is  now  inaccessible.  I  propose  to  attack  this  problem  with  a  combined  hardware  and 
software  approach.  Rrst  the  shot-to-shot  reproducibiiity  should  be  stabilized,  e.g.,  by  spatial  filtering 
of  the  laser  beam,  modest  improvements  in  optics,  and  modestly  improved  mechanical  stability;  after 
this  is  done  it  will  become  productive  to  increase  the  digitization  resolution.  Second,  residual  drift  due 
to  temporal  changes  in  the  laser  intensity  pattern  can  be  minimized  by  periodically  updating  the 
stored  background  picture,  in  contrast  to  the  practice  to  date  of  storing  this  picture  once  at  the 
beginning  of  a  data  run.  Finally,  with  inaeased  signal  amplitude  resolution,  several  approaches  to 
extracting  the  interference  signal  will  be  worth  trying.  These  would  be  primarily  fourier  transform 
methods,  in  which,  perhaps  after  preliminary  processing  to  subtract  the  edge  shadow  signal 
components,  we  would  look  for  the  comb  of  spatial  frequency  peaks  in  the  fourier  spectrum  that 
correspond  to  the  target  diameter. 
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Appendix  A:  The  HUYGENS  Program 


hiqfgsiui.e:  huygans  [axgl  arg2] 

calcttlata  aaiplltudas  uaing  Buygana'  Pxlncipla  as  modifiad  by  Kischoff 

CBLCOXATBS  TBB  RKLXTIVX  OGMPLKZ  AMPLITODB  SXBN  AT  BACH  CBLL  OT  A 
ZJNBAR  CCD  THAT  ZMTIBCIPTS  A  COLLXMATBD  XASXR  BXMf  IHTKBUPTXD  BT  A 
mBB-LIXX  OPAQDI  OBJBCT.  THB  COBBXSPOSIDIMG  SIGNAL  IS  TBB  SOGARX 
MAGMITaOX  or  THB  OOTPGT  OP  THIS  PROGRAM:  (anpl.r^2  4-  aagpl.i^2) 

M.  W.  Siagal  -  Robotics  Inatituta  -  Caxnagia  Mallon  ttalvarsity 
May  29,  1986 

Copysi^t,  tha  author  and  nmytom  PPG  Industrias,  panding  daclalon  as  to 
tha  author's  right  to  dasignata  it  as  baing  in  tha  public  doaiain, 
which  is  whara  ha  wants  it  to  ba. 

(1)  no  argunants:  craata  wavafcns  and  distfcns, 

intagrata  from  -LBHS^RADXUS  to  •fLXMS_RADIUS, 
producing  ou^ut  flla  "unobstru"; 

(2)  two  argunants:  raad  back  wavafcns  and  distfcns, 

subtract  out  tha  ragion  batwaan  (and  including) 
tha  two  argunants; 

(3)  [add  latar] :  if  wavafcns  ara  ^^onatric,  taka  this  into  account; 

(4)  [add  latar] :  if  wavafcns  ara  pura  raal  or  pura  imaginary,  taka  this 
into  account. 

Nota:  It  is  usaful  to  dinanslon  arrays  of  (2*n-l) ,  i.a., 
z[-(n-l) : (n-1) ] ,  as  follows: 

<typa>  *x; 

X  ■  (<typa>  *)calloc(2*n-l,  sizaof (<typa>) )  +  a  -  1; 

Givan: 

(1)  a  Ians  of  radius  UNS_RADIt7S,  from  which  anargas 

(2)  a  collinatad  baan  of  gausslan  profila,  axp(- (y_8rc/y00)  ^^2) , 

(3)  timas  a  parabolic  shading  (1  -  (y_src/UHS_RADIUS)^2) 

(4)  to  naka  a  saiooth  transition  to  zaro  intansity,  and  thus  to  aliminata 
structure  from  tha  diffraction  pattam  of  tha  Ians  apartura; 

and 

(5)  a  ced  consisting  of  CBLLS  calls,  CBLL__SIZB  microns  cantar-to-cantar, 

(6)  locatad  such  that  Cantar_Call  »  atoi (argv[l] )  is  on  tha 
axis  of  tha  laxis,  and  thus  on  tha  cantar  lina  of  tha  baam, 

and 

(7)  a  wira  that  obstructs  tha  baam  DISTAMCB  microns  from  tha  plana  of 
tha  ced, 

(8)  whosa  laft  adga  is  at  LaftJBdga  «  atoi(argv[2] )  microns  from 
tha  axis, 

(9)  and  whosa  right  adga  is  at  Rigbt_Edga  »  atoi (argv[3] )  microns  from 
tha  axis; 

(10)  rOR  BXAMPLB:  huygans  116  -50  150 

calculatas  tha  lina  laiaga  rasponsa  of  tha  ced  whos  116th  call  is  on 
tha  Ians  axis,  whan  obstructed  by  a  wira  whosa  laft  adga 
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1«  50  aicxons  to  th«  loft  of  tho  Ions  axis,  and  whoso  right  odgo  is 
150  sdcroxis  to  tho  sight  of  tho  Ions  axis. 

Tho  program  works  as  follows: 

(1)  physical  and  gocnotsical  parasiotors  of  tho  apparatus  aso  DBrniXod; 

(2)  somo  moso  constants  aro  calculated; 

(3)  tho  wawofront  amplitude  is  specified  at  1  adcron  intervals  from  tho 
lens  axis  to  tho  Ions  radius:  wavofen[IilllS^lUU>IUS] ; 

(4)  a  phase,  1/s  distance,  and  (1  -f  cos  (angle))  obliq^ty  product  function 
is  specified  at  1  micron  spatial  incsosmnts  of  tho  offset  between 
source  point  (y_sre)  and  detector  point  (y_cod) ;  since  this  offset 

can  be  between  seso  (y_ssc  «  y_ccd)  and  UDIS_RhDIXJS  -f  ClLTiS  *  CBLL_SIZK 
microns  (from  the  left  edge  of  the  lens  to  the  ri^t  edge  of  the 
ced  when  CenterjCell  «  0,  or  vice  versa,  from  the  right  edge  of  the 
lens  to  the  left  edge  of  the  ced  when  CenterjCell  *  CKLTiS  -  1)  : 
dlstfcn[UEHS_RU}inS  +  CSLLS  *  CBLL_SIZI]  .  ' 

(5)  the  wavefcn[]  and  distfcn[]  are  independent  of  the  specific  location 
of  the  ced,  os  the  location  of  the  ced  or  the  wise,  in  any  application; 
thus  once  the  arithmetic  of  calculating  these  is  done,  the  results  are 
written  onto  files  whose  names  are  hasd->wised  to  "wavefens”  and 
"distfens"  in  the  directory  from  which  huygens  is  executed; 

(6)  if  the  files  "wavefens"  and  "distfens"  exist  in  the  at  rtin  time,  they 
are  read  into  wavefen[]  and  distfcn[],  thus  eliminating  tedious 
se-ealculation . 

CADTION:  obviously,  I  hope,  if  the  DBriMBs  are  changed  and  the  program 
recompiled,  it  is  iaportant  that  the  files  "wavefens"  and 
"distfens"  are  ABSENT  from  the  local  directory  the  first  tism 
the  recompiled  program  is  sun! 

(7)  for  each  ced  cell,  located  at  y_ccd, 

(8)  integrate  in  one  micron  intervals  over  CXXiL__SIZB  isicrons  in  the 
ced  plane 

(9)  and  in  one  adcron  intervals  from  -LENS__BADIUS  to  -fUHSJRADIUS, 

(10)  except  for  the  values  of  y_src  obstructed  by  the  wise. 

Then  print  the  result  in  some  useful  foraiat,  e.g. , 

an  aaiplitude  file  and  a  signal  file. 

*************************************************************************  f 

iinclude  <stdlo.h> 

finclude  <nath.h> 


fdefine  DISTANCE 

167000.0 

/* 

microns,  wire  to  CCD  */ 

idefine  NAVELENGTH 

0.895 

/* 

ndcrons  */ 

fdefine  LENS_RADIUS 

(int)2S00 

/* 

microns  */ 

fdefine  CELL  SIZE 

(int)13 

/* 

microns  */ 

fdefine  SDB  CELL  SIZE 

(int)5 

/* 

adcrons  */ 

Th«  active  width  of 

each  pixel 

fdefine  CELLS 

(lnt)256 

/********************************************************************* 
The  calculation  will  be  done  between  -(CELLS  -  1)  (CELLS  -  1) 

********************************************************************* ! 
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idafina  MU 
Idafina  IMBM 


0.1 

30.0 


/*  nuaarical  apartuxa  */ 

/*  dag  lasar  divaxganea  */ 


■aln  (axge,  azgv) 

iat  axgc; 

char  *axgr[]; 

{ 

OOMPm  *wava£ea,  *dlst£cn,  *aap; 
/****************«************************************** 

Meha  tYpada£  o£  COMPZXX  should  ba  shxuct  {doubla  x,  i) 
******************«**********«******************«*******/ 

Int  scratch, 

Lana_Radlua  «  (lot)  LBMSJtADZTJS, 

Calls  >  (Int)  CBLLS, 

Call_Slsa  «  (Int)  CBLL^SXZK, 

La£tjBdga  >  atol(axgv[l]) , 

RlghtJBdga  »  atoi (axgv[2] ) , 

y_alsau  «  Lans_Radlu8  -f*  Calls  *  Call_Siza; 
char  *calloc(), 

out_naBa[32] ; 

riU  *£opan(),  *out_£p; 

1£  (axgc  !>  1  M  axgc  !■  3) 

{ 

£pxlnt£ (stdaxx, 

"Usaga :  huygana  [La£t  Bdga  Rlght_Kdga] \n” ) ; 
aaitO; 

) 

scratch  ■  (axgc  «  1)  ?  1  :  0; 

wava£cn  •  (CONPUaC  *)calloc(2*Xians_JUdius  -  1,  sizao£  (CQNPLSX) )  -f 

Xians_Radlua  -  1; 

gat_wa'«a£cns  (wava£cn) ; 

dist£cn  «  (CQMPUEX  *)calloc(2*y_alam8  -  1,  sizao£  (COMPLSX) ) 

yjalsns  -  1; 

gat_diat£cns  (dist£cn) ; 

amp  «  (COMPLEX  *)calloc(2*Calls  -  1,  sizao£ (COMPLEX) )  * 

Calls  -  1; 

calcalata_aBp(waTa£cn,  dist£cn,  La£t_Edgo,  Right_Bdga,  sap,  scratch); 
i£  (scratch)  sprint£  (out_naBia,  "unobstru”); 

alsa  sprint £(out_naina,  "L%<St.%d.aiip'',  La£t_Bdga,  RightJEdga) ; 

out__£p  ■  £opan(out_naaia,  "w"); 

£writa(aap  -  Calls  +  1,  8izao£ (COMPLEX) ,  2*Call8  -  1,  out_£p) ; 

£cl08a (out_£p) ; 

) 

gatjiravo£cn8  (wava£cn) 

COIOLEX  *wava£cn; 

( 

PILE  *wava_£ila,  *£opan(); 

i£  ( (wava  £lla  «  £opan("wava£cns'’,  "r”))  NULL) 

( 

doubla  lans_radlu8  ■  (doubla)  LENS  RADIUS, 
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fwbm  a  (doobl*)  niHM, 

na  >  (doubla)  lA, 

pi,  yOO,  yyO,  yyl,  yy2; 

Int  LwasJRadittS  >  (i&t)  LKHSJADITJS, 

y_8rc; 

pi  a  4.0  *  atan(l.O); 

yOO  a  l«ns__raditts*tan(0.5*£wba*  (pi/180.0)  )*sqrt  (-log (0.5) ) /tan (asln(na) ) 

for  (y_axc  a  -z«aa  Radiua  -f  1;  y_sxe  <  Lana  Radius;  -Hy^src) 

{ 

yyO  a  (doubla)y_ssc; 

yyl  a  yyO/yOO; 

yy2  a  yyO/lans^radius; 

(waTafcn  +  y_8tc)->t  a  axp(-yyl  *  yyl)  *  (1.0  -  yy2  *  yy2) ; 

(wavafen  y  arc)->l  a  o.O; 

> 

wava^fila  a  fopanC’wavafens",  ”w“); 

twsita  (wava£cn-X>ana  Radius^f  1 ,  sixaof  (COMPLRX) ,  2*Lans  Radius-1 ,  wava  f ila) ; 

> 

alsa 

( 

int  Itans^Radius  a  (int)LEHSJRRDIUS; 

fxaad (wavaf en-Lans  Radius4-1, sixaof (CQMPXXX) ,2*Lana  Radiu8-l,wava_£ila) ; 

) 

felosa (vara  fila) ; 

I 

gat_dist£cns (distfen) 

CXRffUEX  *dist£cn; 

{ 

rZLB  *dist_£ila,  *£opan(); 

if  ((dist  fila  a  fopanC'distfcns”,  ”r”))  as  mdll) 

( 

int  Lans^Radius  a  (int)  LBMS_RRDZUS, 

Call?  a  (int)  CBLZJ, 

Call_Sixa  a  (int)  CSLI._SIZX, 

y_alams  a  Lans^Radius  Calls  *  Call_Sixa, 

y  ssc  cod; 

doubla  <Ustanca  -  (doubla)  DISTRMCI, 

wavalangtb  a  (doubla)  WRVXLBISTB, 

pi,  kOO,  pathlan,  dist,  dsq,  cosphi,  sinphi,  kixchf; 


pi  a  4.0  *  Stan (1.0); 
kOO  a  2.0  *  pi  /  wavalangth; 
dsq  a  distanca  *  distanca; 

for  (y  sxc  ced  »  0;  y_8rc_ccd  <  y_alans;  +4y  src_ccd) 
{  “  ” 

dist  a  (doubla)  y__src_ccd; 

patblan  a  sqrt  (dsq  +  dist  *  dist)  ; 

cosphi  a  cos (kOO  *  patblan) ; 

sinphi  a  sin (kOO  *  pathlan) ; 

kixchf  a  (1.0  +  distanca/pathlan) /patblan; 
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(distfon  y_scc__ccd) ->r  «  (dlstfen  -  y_arc_ccd) ->c  *  coaphi  *  klrchf; 
(dlatfcn  -f  y__arcjeed) ->1  «  (dlatfen  -  y  arc_jeed)->i  ■  ainphi  *  kirchf; 
) 

diat_£il«  >  fopanC^diatfcna",  "w"); 

fwrltai  (diat£ea  '  y_«laiaa  4-  1,  alzao£  (COMPLEX) ,  2*y  alana-l,  diat_£ila); 

) 

•laa 

( 

int  Lana^IUdiua  «  (int)  LENSJRADIUS, 

Calla  »  (Int)  CELLS, 

Call_Siz«  -  (int)  CELL_SZZE, 

yjalaoa  «  Lana__Iladiua  Calla  *  C«ll_Sla«; 

£z«ad(dlat£en  -  y_«laima  4-  1,  aiaao£  (COMPLEX) ,  2*y  alaaw  -  1,  diat_£il«)  ; 

) 

£cloa« (diat_£ll«) ; 

) 

cnleulatai_ai!^(wava£cn,  diat£en,  La£tJBdga,  Ri^t_Edga,  aap,  acratch) 
COMPUDC  *wa'v«£cn,  *diat£cn,  *aiip; 

Int  La£t_Xdga,  Right  Edga,  acratch; 

{ 

int  Calla  a  (int)  CELLS, 

Call_Siaa  >  (int)  CELL_SIZE, 

SubjCall_Siza  a  (int)  SCBJOCLL  SIZE, 

Sab_La£t  a  -(Sub_Call_Siaa/2) , 

SubJRight  a  Sub_Call__Sisa/2  +  1, 

LaaaJUdiua  a  (intT  LEHS_RADins, 

call,  aubcall,  y_ccd,  y_arc; 

COMPLEX  aaf>l,  *d_y_ccd;  ~  *” 

1£  (acratch) 

( 

La£t__Edga  a  -LanaJRadiua  4-1; 

Rig^t_Edga  a  Lana  Radiua  -  1; 

} 

alaa 

( 

FIU  *unobatru__£p,  *£opan(); 
unobatru__£p  a  £opan(”unobatru",  '*r'*); 

£raad(aBg>  -  Calla  4-  1,  aizao£  (C^fPLEX) ,  2*Calla  -  1,  unobatru  £p) ; 

} 

£or  (call  a  -calla  4-  1;  call  <  Calla;  4-l-call) 

{ 

aapl.r  a  o.O; 
aapl.i  a  0.0; 
y_jccd  a  call  *  Call_Siza; 

£or  (aubcall  a  Sub_La£t;  aubcall  <  Sub  Right;  4-faubcall) 

I 

d_y_ccd  a  diat£cn  -  y_ccd  -  aubcall; 

£or  (y  arc  a  La£t  Edga;  y_8rc  <a  Right  Edga;  4-fy_arc) 

aapl.r  +a  (d_y_ccd  +  y_arc)->r  *  (wava£cn  +  y_arc) ->r  - 
(d_y_ccd  4-  y_arc)~>i  *  (wava£cn  +  y_arc)->i; 
aapl.i  +a  (d_y_ccd  4-  y_arc)->r  *  (wava£cn  +  y_arc)->i  4- 
{d_y_ccd  +  y_arc) ->i  *  (wava£cn  +  y_arc) ->r; 
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} 

) 

if  (scratch) 

( 

lamp  call)->c  ■  anpl.r; 
(•flip  -f  call)>>i  -  aat>l.i; 

> 

alsa 

( 

(asp  +  call)->r  -■  ain>l.r; 
(aa^  +  call}'->i  -■  aapl.i; 
> 

} 

I 
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Appendix  B:  The  INCLUDE  FWe  of  Constants  and  Definitions 


iincluda  <«tdlo.h> 
ilaelud*  <Bath.h> 


idafliM 

BXG6 

1.0*9 

#d»£iM 

sinr 

1.0*-9 

idAfln* 

PI 

3.141592654 

idaflna 

CCDS 

4 

idafln* 

rUT.T.g 

256 

idafixM 

AMPCXLIiS 

511 

/* 

£or  making  th*  amsk  £raia  huyg*ns  ( 

#d*£iM 

AOC_IUIJ«GB 

63 

/* 

63  dark,  0  saturatad  light  */ 

#d*£lxM 

MIN_ABXX 

1.0 

/* 

sqoara-ccds  */ 

#d«£in« 

DIXJBSR 

10.0 

/* 

ccd  */ 

idafln* 

CBI_BRR 

2.0 

/* 

ccd  */ 

idaflxM 

MXSKLBH 

64 

/* 

convolution  mask  langth  */ 

idafliM 

MXSKZKR 

2 

/* 

Biask  o££s*t  */ 

idafln* 

HOirjSDTSY 

16 

/* 

willingnass  to  risk  */ 

/* 

intgr  math  ov*r£low  */ 

fdafixM 

imiGBT 

0.001 

/* 

ralativ*  waight  o£  cantar  */ 

/* 

vs.  diamatar  arrors  in  */ 

- 

/* 

tka  corraction  algorithm  */ 

fdAflXM 

aU._TABLB 

"cal" 

/* 

root  nama  o£  calibration  £ilas  */ 

IdafliM 

OOnrABLB 

"c*n" 

/• 

root  nama  o£  cantaring  £ilas 

*/ 

idafliM 

tmoisTKD  rzLS 

’•nif" 

/* 

"no  wira"  */ 

id*£ln« 

MOj:ASER~rZlX 

"nl" 

/* 

"no  lasar"  */ 

idaflxi* 

Mnikjrzu 

"mask" 

/* 

naad  I  say  mora?  */ 

idaflxM 

I*Y 

••/dsv/ttyOO” 

fdafln* 

BAXJD 

B9600 

idaflM 

XNQ 

5 

/* 

CELLS  data  points  sant  on  racaipt 

fd«£ln« 

SPACING 

13.0 

/*  microns  */ 

#d*£lu 

DISTANCB 

167000. 

0 

/*  microns  */ 

id«£in« 

GK_PBR_CK3 

2.53 

#d*£ln« 

PACKING__rRAC 

0.91 

idmtlam 

GItAPH_NONX 

0 

/* 

thasa  nin*  ar*  anum-s  in  tha 

*/ 

tdmtinm 

6RAPH~S(MB 

1 

/* 

main  program, 

*/ 

idmiiam 

6I(APH_AZJi 

2 

/* 

but  I  had  di££iculty  passing 

*/ 

/* 

them  to  subroutinas  as  such 

*! 

idmtlnm 

PAUSX  NONE 

0 

/* 

and  rathar  that  Tight  it  now 

*! 

#da£ln* 

PAUSE  SOME 

1 

/* 

I  pass  tham  as  intagars 

*! 

ida£ln« 

PAUSEJOJi 

2 

#da£ln« 

CAZiZB_I.IVE 

0 

#da£ln« 

CAZJB  AUTO 

1 

#da£ln« 

CAZ.IBrMANU 

2 

struct 

DC  <£10St 

d,  c; )  ; 

struct 

POINT  (£10St 

*/  y;); 

striict 

LINE  {£lost 

m,  b;); 
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struct  PLANS  (float  cx,  cy,  c*;>;  /*  *  ■  cac*x  +  cy*y  +  c*  */ 

struct  PTPAIR  (struct  POINT  corr,  bmss;};  /*  corrsct  and  aaasurad  */ 

struct  DPP  (float  dlst;  /*  (corr  -  aaas)^2  */ 

struct  POINT  corr,  mss;); 

typadaf  struct  (double  r,  i;}  COMPLEX; 
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Appendix  C:  The  DIACENS  Diameter-Center  Aigorithm  Program 


/* 

dlac*na() 

glvvn  (char)  cdata,  xatuma  (dc)  dlaiaatar-cantar  pairs 

*/ 


/****«**********«*********************«*********************************** 

*  N.  W.  Slagal  -  Robotics  Institute  -  Casnagla  Mellon  University  * 

*  * 

*  Copyright  1987  and  earlier  developsiant  versions,  the  author  and  * 

*  possibly  PPG  Industries,  pending  decision  as  to  the  author's  right  * 

*  to  deslgxiate  It  as  being  In  the  public  donaln,  which  Is  where  he  * 

*  wants  it  to  be.  * 


tlnclude  "yardage. h" 

dlacens  (edata,  dc,  graph_aode,  patise_mode,  eallb_mode) 
unsigned  char  *cdata; 

struct  DC  *dc; 

Int  graph_8K>de,  pause  aoda,  eallb  xaode; 

{ 

register  Int  leng,  conv,  *lnte,  *iBask,  cell,  convmln; 
char  caption [80 ] ; 

Int  ced,  rcell,  rmln,  Icell,  Imln, 

IMTKMS  [CELLS], 

COMV  [CELLS] , 

VMOC  [CELLS] ; 

double  sedge,  ledge,  centr,  apeak (); 

struct  DC  precalO; 

static  double  radix; 

static  int  nevercalled  •  1, 

MhSK  [AMPCELLS] ; 

Static  float  NOJEIRE [CCDS] [CELLS] , 

NO_IASE[CCDS] [CELLS] , 

nramis  [cells]; 

If  (nevercalled) 

{ 

riU  •fp_in; 
float  ftenp; 

nevercalled  «  0; 

radix  a  hOWjSUTSY  *  /*  H0W_G0TSY:  willingness  to  risk  overflow  */ 

pow(2.0,  0.5  *  8.0  *  slzeof(lnt))  / 

(ADC_RAMGE  *  sqrt ( (double) (CELLS  -  MhSKLEN)  *  MESKLEN) ) ; 

/*  get  and  optionally  show  the  mask  */ 

fp_in  ■  fopen(MESK_PILE,  "r") ; 

for  (coll  «  0;  cell  <  AMPCELLS;  -H-cell) 

{ 

fscanf (fp_in,  "%f",  fiftenp) ; 

NASX[cell]  «  fteap  *  radix  ■¥  0.5; 

) 

fclose(fp_ln) ; 
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if  (graph  aoda  GRAPH_AIiL) 

{ 

graph («MhSK[AIIPCBLX.S/2],  'i',  AMPCSLLS/2,  1,  ''mask  (right  half) ") 

await_cr (pausa  noda) ; 

> 

/*  gat  and  optionally  show  tha  baekgrouxid  (no_lasar)  fila  */ 
fpJLn  >  fopan(MO_LhSKR_riLS,  "r"); 
for  (ced  «  0;  ccd  <  CCDS;  -Kfccd) 

for  (call  w  0;  call  <  CBU<S;  •(-Kwll) 

faeanf  (fpJLn,  '*%f " ,  SMO^LASB [ccd]  [call] )  ; 
fclosa (fp_in) ; 

if  (graph_aoda  »  (3tAPa_ALIi} 

for  (ccd  w  0;  ccd  <  CCDS;  -H^cd) 

( 

sprintf (caption,  "ccd  %ld:  no_^lasar  fila",  ccd); 
graph (MO_LASB [ccd],  'f',  CELLS,  1,  caption); 

await_cr  (pau8a_inoda) ; 

> 

/*  gat  and  optionally  show  tha  unobstructad  (nojwira)  fila  */ 
fp_in  ■  fopan(UMOBSTBD_riLB,  "r“); 
for  (ccd  >  0;  ccd  <  CCDS;  -H-ccd) 

for  (call  a  0;  call  <  CELLS;  >-K:all) 

fscanf (fp_in,  "%f",  SNOJTCRS [ccd]  [call] ) ; 
fclosa (fp_in) ; 

if  (graphjBoda  6RhPH_ALL) 

for  (ccd  w  0;  ccd  <  CCDS;  'f+ccd) 

( 

sprintf (caption,  "ccd  %ld:  no_wira  fila",  ccd); 
graph (HOJEIRE [ccd],  'f',  CELLS,  1,  caption); 

await_er (pausa  aoda); 

) 

) 


/* 

EAENIEG:  doas  not  includa  any  NO_LASER  corractions, 
bacausa  I  am  not  raally  sura  how  it  should  ba  dona, 
or  awan  if  it  should  ba  dona  at  all! 

*/ 

for  (ccd  a  0;  ccd  <  CCDS;  -M-ccd,  -Hdc,  cdata  +m  CELLS) 

( 

/*  SBiooth,  convert  to  floats,  nosaaliza  tha  input  data  (chars)  */ 
f8aooth( cdata,  rMTENS,  graph_iDOda,  pausa_noda,  ccd); 

/*  correct  for  bac)cgrouzul  */ 

subnozm(rNTENS,  NO_inRB[ccd] ,  graph_aode,  pausa_aK>da,  ccd)  ; 

/*  convert  to  integer  arithmatic  */ 

for  (cell«0;  calKCELLS;  -H-call)  IIlTEMS[call]«niTENS[call]*radix+0.5; 

/*  taka  convolutions  to  find  tha  right  edge  */ 
for  (call  «  0;  call  <  CELLS  -  MhSKLEM;  -H-call) 

{ 

conv  »  0; 
lang  -  MESKLEM; 

aask  >  &MASK[AMPCELLS/2]  -I-  MESKZER; 
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lnt«  >  ZNTENS  +  call; 

whlla  (  laxig —  )  conv  +■  *inta++  *  ^nask-H-; 

OQNVtcall]  •  conv; 

> 

for  (call  ■  CELT'S  -  MASKLEM;  call  <  CELLS;  -H-call)  CONV[call]  «  conv; 
convfldA  «  CONV[01 ; 

for  (call  >  1;  call  <  CELLS  -  KASXLEM;  -H-call) 

If  (COIlV[call]  <  convodn)  convxnin  »  COMV[xmin  «  call]  ; 
call  >■  (reall  *  rmin)  -f  1; 

whlla  (COMV[rcall]  <  CONV[call-H-] )  -H-rcall; 
radga  «  rcall  -f  zpaa)c (&COMV[ reall]) ; 

/*  taka  convolutions  to  find  tha  laft  adga  */ 
for  (call  >■  MASKLEM;  call  <  CELLS;  -H-cell) 

( 

conv  a  0; 

lang  ■  MASKLEN;  /*call<KASKLEK  ?  call  :  MASXLEN;*/ 
mask  a  CMASK[AMPCELLS/2]  -t-  NASKZER; 
inta  a  INTEHS  ■¥  call; 

la  (  lang —  )  conv  -fa  *lnta —  *  *inaak-H-; 

VMOCCcall]  a  conv; 

I 

conv  a  VMOC(MASKLEH] ; 

for  (call  a  0;  call  <  MASKLEM;  -H-call)  VMOC[call]  a  conv; 
convBiin  a  VHOC[0]; 

for  (call  a  maSKLEH;  call  <  CELLS;  -H-call) 

if  (VMOC[call]  <  convmin)  convnin  a  VNOC[lmin  a  call]; 
call  a  (lean  a  lain)  -  1;  whila  (VMOC[lcall]  <  VNOC[call— ])  — Icall; 
ladga  a  icall  -f  zpask(SVNOC [Icall]); 

/*  do  soma  spacial  stuff  if  closa  to  a  ced  and  */ 

if  (  min  <a  maskleK  )  /*  than  tha  LEFT  adga  is  garbaga  */ 

( 

dc->c  a  0.5  *  (rmin  -f  zpaak(fiCONV[nnin] )  -f  lsd.n  -f  xpaak(SVNOC[laiin] ) ) 
dc->d  a  2.0  *  (radga  -  dc->c) ; 

*dc  a  pracal(*dc); 

> 

alsa  if  (  Imin  >  CELIiS  -  NASKLEN)  /*  than  tha  RIGHT  adga  is  garbaga  */ 

( 

dc->c  a  0.5  *  (rain  +  zpaak(fiCOMV[xmin] )  -f  Loin  +  zpaak(&VNOC[lmin] ) ) 
dc->d  a  2 . 0  *  (dc->c  -  ladga) ; 

*dc  a  pracal (*dc) ; 

) 

alsa  /*  this  is  tha  normal  casa  */ 

{ 

dc->c  a  0.5  *  (radga  -f  l-iga)  ; 
dc->d  a  radga  -  ladga; 

*dc  a  pracal (*dc); 

) 

cantr  a  dc->c; 

if  (  graph  moda  |  |  calib  moda  as  CALIB_MAMU  ) 

( 

sprintf (caption,  "ced  %ld:  diaaatar  %7.3f,  cantar  %7.3f", 
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ccd,  dc-Xi,  dc->c)  ; 

gxaph(  VHOC,  'i',  CBtLS,  1, 

gxaph(  CONV,  'i',  CELLS,  0,  caption); 

graphjBarlc  ( (float)  ladga,  -1,  '  >' ) ; 
grapljMBark  ( (float)  cantr,  -1 ,  '  | ; 
graph_maxk ( (float) radga,  -1 ,  ' <' ) ; 
await_er  (pauaa_]aoda) ; 
if  (  eallbjnodia  CALIB  MAUD  ) 

{ 

graph(INTENS,  '1',  CXLLS,  0,  caption) ; 

aanu  cal(*dc,  ccd); 

> 

> 

)/*  And  of  ccd  counting  loop  *! 
ratuxn; 

> 


fsaooth(cdata,  intana,  gxaph_aioda,  pausa_i&oda,  ccd) 


/*  sa»oth  character  data  to  float,  and  perhaps  graph  it  */ 

/**********************************************************************/ 
unsigned  char  *cdata; 
float  *intens; 

int  graph_mode,  pause  mode,  ccd; 

{ 

int  cell ; 


*(intens  +  0)  *  *(cdata  +  0); 

*(intens  +  1)  *  (*(cdata  +  0)  + 

2.0  *  *(cdata  +  1)  + 

*(cdata  <f  2))/4.0; 

for  (cell  ■  2;  cell  <  CELLS  -  2;  -H-cell) 

*(intens  +  cell)  ■  (*(cdata  +  cell  -  2)  + 

2.0  *  *(cdata  ■¥  cell  >  1)  -1- 

3.0  *  *(cdata  +  cell  )  + 

2.0  *  *(cdata  +  cell  +  1)  + 

*(cdata  +  cell  +  2))/9.0; 
*(intens  -f  CELLS  -  2)  »  (*(cdata  4-  CELLS  -  3)  -t- 

2.0  *  *(cdata  4-  CELLS  -2)4- 

*(cdata  4-  CELLS  -  l))/4.0; 
*(intens  4-  CELLS  -  1)  »  *(cdata  4-  CELLS  -  1); 


if  (graph_mode  sa  GRAPH  ALL) 

( 

char  caption [80 ] ; 

sprintf (caption,  "ccd  %ld:  smoothed  data” ,  ccd) ; 
graph (intens,  'f',  CELLS,  1,  caption) ; 

await  cr(pause_mode) ; 

) 

return; 

) 


subnozm (intens,  nowira,  graph_moda,  pause_mode,  ccd) 
^********************************^**11^*****^ ******* **11^******** **********  ^ 

/*  subtract  from  a  background  pattern,  and  perhaps  graph  it  */ 
^************************* ****************************** *************** ! 
float  *intens,  *nowire; 
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Int  graph  aoda,  pauaa_noda,  eed; 

{ 

raglatar  int  call; 

call  «  CELLS; 
whlla  (  call —  ) 

{ 

*lntans  •  *nowira  -  *lntans; 

-H-intans;  -H-nowira; 

} 

iataaa  CELLS; 
if  (graphjaoda  GRAPH  ALL) 
char  caption [80]; 

aprlntf  (caption,  "ccd  %ld:  background  axibtractad  and  noxmalizad" ,  ccd)  ; 
graph (intaxia,  'f',  CELLS,  1,  caption); 

await_cr (pauaa  mode) ; 

> 

ratum; 

} 

*/ 

zpaakO  :  intarpolation  for  aaxiimiiB  or  miniaam  in  an  array 
dariv() :  darivativaa 
fib() :  fibonaccl  nunbars 

*/ 

doubla  zpaak(y) 
int  *y; 

/*******«*********«****************************************************** 
Glvan  *(y  -  1)  <  •y  >  *(y  +  1)  or  *(y  -  1)  >  *y  <  *{y  +  1) , 
rpaak(y)  ratuma  tha  valua  of  z  (-1.0  <  z  <  1.0)  at  which  tha 
maTlwBim  or  miniiwm  "actually”  occurs,  basad  on  an  azpansion 
fitting  tha  ragion  undar  tha  curva  from  *  (y  -  DY)  to  *  (y  DY) , 
i.a.,  it  will  consider  deriwativas  of  y  up  to  d(2*DY)y/dz(2*DY) 
and  powers  of  z  up  to  z^ (2*DY) . 

Tha  procedure  (Newton's  z(n+l)  >  z(n)  -  y' (n)/y'' (n) )  is  iterated  until 
suoqassiva  values  of  z  differ  by  lass  than  DX. 

************************************************************************* ^ 

{ 

fdafina  OX  0.001 
idefina  MI  6 

doubla  yl,  y2,  y3,  y4,  y5,  y6,  a,  b,  z,  zl,  z2,  z3,  z4,  dariv(); 
int  itars; 


yl  ■  dariv(y. 

1, 

1.0) 

y2  >  deriv(y. 

2, 

1.0) 

y3  a  dariv(y, 

3, 

1.0) 

y4  ■  dariv(y, 

4, 

1.0) 

y5  ■  dariv(y, 

5, 

1.0) 

y6  «  dariv(y. 

6, 

1.0) 
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1£  (  yl  —  0.0  II  y2  n  0.0  )  nturn  0.0; 

X  ■  -yl/y2; 

Itxrs  a  MI; 
do  ( 

xl  ■  x; 
x2  a  zl*x; 
x3  >■  x2*x; 
x4  «  x3*x; 

if  ((•  >  yl  y2*xl  +  y3*x2  +  y4*x3  +  y5*x4)  —  0.0)  broak; 
if  ((b  »  y2  +  y3*xl  +  y4*x2  *  y5*x3  +  yfi*x4)  —  0.0)  bxoak; 

X  -«  a/b; 

)  whila  (  faba  (x  -  xl)  >  DX  &£  itars—  ) ; 
ratuzn  x; 

) 

/*  9anaxalized  varsion  of  xpaak() : 

{ 

idafina  DY  3 
idafina  DX  0.001 

doubla  X,  dx,  daziv(),  taylorO,  dar[2*DY  +  1],  tayl,  tay2;  int  i; 
for  (i  a  0;  1  <a  2*OY;  -H-l)  daz[l]  a  daziv(y,  1,  1.0); 

X  a  dx  a  0.0; 

whila  (  (tay2  a  tayloz(daz,  2*Dy-l-l,  2,  x))  ! a  0.0  44 

faba(dx  a  -taylor(daz,  2*OY-l-l,  1,  x)/tay2)  >  OX  )  x  -f*  dx; 

zatuzn  x; 

) 

*/ 

doubla  dariv(y,  n,  dx) 
int  *y; 

int  n; 

doubla  dx; 

/************************************************************************* 
TaJcaa  y,  a  pointar  to  an  alanant  of  an  array  of  doublas,  and 
ratuma  tha  n-th  darivativa  in  tba  vicinity  of  that  alamant. 

Tha  alamants  of  y  ara  saparatad  by  dx,  i . a . ,  tha  horizontal 
dlstanca  batwaan  (y  N)  and  (y  N  -  1}  is  dx. 

if  (  dx<a  0.0  dx»  1.0  )  then  tha  result  is  reported  without 
normalizing  by  dx'^n,  auking  tba  arithmetic  a  little  faster. 

RKQOIBXS:  fib(n,  m) ,  a  Fibonacci  number  generator. 

MkRlIING:  it  is  tha  caller's  responsibility  that  (y  -  (n-fl)/2)  through 
(y  +  (n-l'l)/2)  be  valid  addresses  containing  mambars  of  array  y. 

WARMING:  severe  rounding  errors  may  start  to  appear  around  n  ^  8. 
********************************«***«************************************/ 
{ 

int  *yl,  *yr; 
doubla  dny__dxn; 

int  flb(),  m  ■  0,  sign  »  1,  n_2; 

if  (  n  <  0  )  arr_axit ("dariv:  caimot  taka  derivative  of  order  <  0."); 
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If  (  n  av  0  )  dny_dxn  a  *y; 


if  (  n  %  2  )  /*  odd  a,  oven  nmaber  of  points,  so  syametrizs  */ 

{ 

n_2  a  (n+l)/2; 

dojf_dxa  a  (doubls)  *  (yr  a  y  +  n__2)  -  (doubla)  *  (yl  ay-  n_2)  ; 
whlla  (sign  a  -sign,  -Mib  <  nj2) 

dny_dxn  sign  *  (fib(n,  m  -  1)  -  fib(n,  m))  * 

( (doubla)  *++yl  -  (doubla)  yr) ; 

dnyjdxn  *a  0.5; 

> 

alsa  /*  avan  n,  odd  nuabas  of  points,  is  syssaatric  */ 

{ 

n_2  a  n/2; 

dny_dxn  a  (doubla)  *  (yl  ay-  n_j2)  +  (doubla)  *  (yr  a  y  +  n_2)  ; 
wbila  (sign  a  -sign,  -t-Hs  <  n_2) 

dny_dxn  +a  sign  *  fib  (n,  a)  *  ( (doubla)  *-Hyl  +  (doubla)  * — yr)  ; 
dny  dxn  -fa  sign  *  fib(n,  a)  *  (doubla)  *y; 

) 

) 

if  (dx>0  ££dx  !ai.o  )  whila  (  n —  )  dny_dxn  /a  dx; 
ratum  dny__dxn; 

> 


int  fib(n,  m) 
int  n,  a; 

y ************************************************************************** 
Fibonacci  nusbar  ganarator;  on  tha  assuaption  that  tha  first  faw 
▼aluas  ara  naadad  aoat  of tan,  thasa  ara  storad  in  a  static  array. 

**************************************************************************  j 


{ 

idafina  BIGHT  8 

static  int  Fib [BIGHT] [BIGHT]  »  (( 

{ 

{ 

{ 

( 

{ 

( 

( 


1, 

0, 

0, 

0, 

0, 

0, 

0, 

0), 

1, 

1, 

0, 

0, 

0, 

0, 

0, 

0), 

1, 

2, 

1, 

0, 

0, 

0, 

0, 

0), 

1, 

3, 

3, 

1, 

0, 

0, 

0, 

0), 

1, 

4, 

6, 

4, 

1, 

0, 

0, 

0), 

1, 

5, 

10, 

10, 

5, 

1, 

0, 

0), 

1, 

6, 

15, 

20, 

15, 

6, 

1, 

0), 

1, 

7, 

21, 

35, 

35, 

21, 

7, 

1}); 

if  (n  <  0  II  a  >  n)  arr  exit ("fib:  called  with  n  <  0  or  a  >  n.") ; 


ratum  (a  a  0  |  |  a 
(n  <  EIGHT 


> 


n)  ?  1 

)  ?  Fib[n]  [a]  : 

fib(n-l,  m-1)  +  fib(n-l,  a) ; 


struct  DC  pracal(dc) 
struct  DC  dc; 

( 

struct  DC  dc_; 

doubla  xl,  x2,  x3,  x4,  yl,  y2,  cO,  cl,  c2,  c3,  c4; 
yl  «  dc.c; 
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y2  ■  yl*yl; 


cO 

m 

-2. 7583*01 

+ 

2.1288*-01  * 

yi 

- 

8.3575*-04 

*  y2 

cl 

m 

1.5949*00 

- 

8.1579*-03  * 

yi 

+ 

3.2038*-05 

*  y2 

c2 

m 

0.0 

+ 

0.0  * 

yi 

+ 

0.0 

*  y2 

c3 

m 

0.0 

+ 

0.0  * 

yi 

+ 

0.0 

*  y2 

C4 

* 

0.0 

+ 

0.0  * 

yi 

+ 

0.0 

*  y2 

xl 

m 

dc.d; 

s2 

m 

xl*xl; 

z3 

m 

xl*z2; 

xA 

m 

xl*z3; 

dc 

.d  a  CO  +  Cl 

* 

xl  c2  *  X2 

+  e3 

*  x3  +  c4  ' 

*  x4; 

.c  >B  dc.c; 

xatuzn  dc  ; 

> 
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(19  cont'd  ) 

merit  are  measured  in  a  tdata  run  of  1000  laser  shots.  Location  and  diameter  measurement 
standard  deviations  are  0.13  pixel-widths  (1.69  pm)  and  0.36  pixel-widths  (4.68  pm) 
respectively,  implying  that  sub-micron  precisions  are  obtainable  by  averaging  a  few  dozen 
measurements.  Sources  of  drift  and  scatter,  and  their  reduction,  are  discussed  in  the 
context  of  designing  and  building  an  improved  second  generation  instrument. 
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